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ABSTRACT 

We present Spitzer IRS spectra (R ~ 600, 10 - 38 jum) of 38 positions in the 
Galactic Center (GC), all at the same Galactic longitude and spanning ±0.3° 
in latitude. Our positions include the Arches Cluster, the Arched Filaments, 
regions near the Quintuplet Cluster, the "Bubble" lying along the same line-of- 
sight as the molecular cloud GO. 11— 0.11, and the diffuse interstellar gas along 
the line-of-sight at higher Galactic latitudes. From measurements of the [O IV], 
[Ne II], [Ne III], [Si II], [S III], [S IV], [Fe II], [Fe III], and H 2 S(0), S(l), and 
S(2) lines we determine the gas excitation and ionic abundance ratios. The 
Ne/H and S/H abundance ratios are ~ 1.6 times that of the Orion Nebula. The 
main source of excitation is photoionization, with the Arches Cluster ionizing 
the Arched Filaments and the Quintuplet Cluster ionizing the gas nearby and 
at lower Galactic latitudes including the far side of the Bubble. In addition, 
strong shocks ionize gas to +3 and destroy dust grains, releasing iron into the 
gas phase (Fe/H~ 1.3 x 10~ 6 in the Arched Filaments and Fe/H~ 8.8 x 10~ 6 in 
the Bubble). The shock effects are particularly noticeable in the center of the 
Bubble, but +3 is present in all positions. We suggest that the shocks are due to 
the winds from the Quintuplet Cluster Wolf-Rayet stars. On the other hand, the 
H2 line ratios can be explained with mult i- component models of warm molecular 
gas in photodissociation regions without the need for H 2 production in shocks. 

Subject headings: Galaxy: center — infrared: ISM — ISM: abundances — ISM: 
bubbles — H II regions 
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1. Introduction 

The Milky Way's Galactic Center (GC) includes a quiescent, ~4 x 10 6 M Q black hole 
(Ghez et al. 2005; Eisenhauer et al. 2005), three clusters of young massive stars, and massive 
molecular clouds (see the reviews in Schodel et al. 2006). Figure la shows a radio image of 
the region from Yusef-Zadeh & Morris (1987a), where the labels refer to the objects discussed 
in this paper. The Radio Arc (Yusef-Zadeh et al. 1984) consists of non-thermally emitting 
linear filaments perpendicular to the Galactic plane; the Arched Filaments and Sickle, on the 
other hand, are thermal emitters, as is seen from their production of radio recombination 
lines (e.g., Yusef-Zadeh & Morris 1987b; Morris & Yusef-Zadeh 1989; Lang et al. 1997, 
2001). Figure lb shows a Band E (21 /im) image from the Midcourse Space Experiment, 
MSX (Egan et al. 1998; Price et al. 2001). At 21 fim, most of the emission comes from 
warm dust. 

Sgr A West is an H II region at the very center of the Galaxy, to which we assume 
a distance of 8 kpc. It contains a cluster of massive stars and the black hole, which is 
coincident with radio source Sgr A*. The other two clusters of massive stars, the Arches 
Cluster and the Quintuplet Cluster, are located about 25 pc away in the plane of the sky. 
The Arches Cluster, whose O and Wolf-Rayet (WR) stars were first discussed by Nagata et 
al. (1995) and Cotera et al. (1996), contains ~ 10 4 M (Stolte et al. 2002), emits - 10 5L6 
ionizing photons s" 1 (Figer et al. 2002), and is suggested to have an age of ~ 2.5 Myr (Figer 
et al. 2002). The Quintuplet Cluster, named after the five luminous near-infrared (NIR) 
objects that may be dusty carbon WR stars (see the discussion by Moneti et al. 2001), is a 
little smaller (10 50 ' 9 s _1 ionizing photons, Figer et al. 1999b) than the Arches Cluster and 
probably older (Figer et al. 1999a, 1999b). 

Given that it takes < 1 Myr for a cluster to orbit the GC, it is probable that either 
or both the Arches and Quintuplet Clusters are now separated from their natal molecular 
clouds. Any interactions with the GC's massive molecular clouds come from the long range 
effects of the clusters' high ionizing luminosities and the strong winds from their massive O 
and WR stars. With discovery of the hot stars of the Arches and Quintuplet Clusters, it 
was immediately suggested that these clusters provide the ionizing photons needed to power 
the Arched Filaments and Sickle, respectively (e.g., Cotera et al. 1996; Figer et al. 1995; 
Timmermann et al. 1996). 

Since the thermal ionized gas was discovered before the star clusters, there were orig- 
inally questions regarding the source of excitation of the gas (e.g., Erickson et al. 1991). 
With the discovery of the two clusters, proof was sought that these are the exciting clusters 
by measuring the excitation variations (i.e., the ionization structure) as a function of dis- 
tance from the clusters. To do this, far infrared (FIR) lines were observed using the Kuiper 
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Airborne Observatory (KAO) (Colgan et al. 1996; Simpson et al. 1997) and the Infrared 
Space Observatory (ISO) (Rodriguez- Fernandez et al. 2001b; Cotera et al. 2005). For both 
the Arched Filaments (Colgan et al. 1996; Cotera et al. 2005) and the Sickle (Simpson et 
al. 1997; Rodriguez- Fernandez et al. 2001b), the excitation decreases with distance from 
the Arches Cluster and Quintuplet Cluster, respectively, as expected for photoionization. 

However, the radial velocity structure of both gas and stars in this region is complex, and 
there could be substantial line-of-sight distances between the star clusters and the ionized gas. 
The Arched Filaments have radial velocities ranging from ~ —25 km s _1 at the northern end 
to rsj —40 km s -1 closest to Sgr A with occasional Vlsr ~ km s _1 at intermediate locations 
(Lang et al. 2001). At many locations the observed lines show more than one velocity 
component. Because negative radial velocities are contrary to the direction of Galactic 
rotation, Lang et al. (2001) suggested that the gas is moving on an elongated orbit in the 
Galaxy's barred spiral potential. The gas of the Sickle has an even larger range of velocities 
with Vlsr ranging from ~ —40 km s _1 to ~ 140 km s _1 (Yusef-Zadeh et al. 1997; Lang 
et al. 1997). Unlike the Arched Filaments, the average Vlsr of the Sickle is positive, ~ 35 
km s _1 (Lang et al. 1997). In both the Arched Filament and the Sickle regions, there exist 
molecular clouds with velocities similar to those of the ionized gas (e.g., Serabyn & Glisten 
1987, 1991; Lang et al. 2002). The ionized gas is most likely the photoionized edges of these 
clouds (Lang et al. 2001). 

Even though it seems likely that the thermal gas in Fig. la is the photoionized edges 
of molecular clouds, the morphology seen in Fig. 1 suggests that the cloud structures are 
not randomly placed in the GC. In particular, the systems of semi-circular arcs suggest that 
the morphology has been influenced by stellar winds or supernova explosions. The most 
notable pair of arcs has become known as the "Bubble" (Levine et al. 1999) or the "Radio 
Arc Bubble" (Rodriguez- Fernandez et al. 2001b). This is the circular open structure south 
of the Quintuplet Cluster whose rim is faint at radio frequencies (Fig. la) but prominent in 
the 21 fim continuum produced by warm dust (Fig. lb, Egan et al. 1998; Price et al. 2001). 

At the same position on the plane of the sky as the Bubble, there is a dense molecular 
cloud called either GO. 11-0.11 (Tsuboi et al. 1997) or GO. 13-0.13 (Oka et al. 2001). The 
column densities through the cloud are so high that the cloud is optically thick even at mid- 
infrared (MIR) wavelengths (A v ~ 640 - 740, Handa et al. 2006). In fact, the GO. 11 -0.11 
cloud is so optically thick that Handa et al. (2006) suggest that the dark area of the center 
of the Bubble is the shadow of the cloud. In this scenario the Bubble rim would be the 
illuminated exterior of the GO. 11— 0.11 cloud. We note, though, that there is no hint of such 
an optically thick cloud in the Spitzer Space Telescope IRAC images of the region (Stolovy 
et al. 2006, 2007), at which wavelengths (3-8 fim) the cloud would be even more optically 
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thick. We question whether the Bubble is in fact related to the GO. 11— 0.11 cloud. We 
suspect the Bubble is a low density, interstellar bubble of hot shocked gas, lying between the 
Earth and the optically thick GO. 11— 0.11 cloud in the GC. 

In this paper we present Spitzer Space Telescope observations (Program 3295, all AORs) 
of MIR spectra of 38 positions along a line approximately perpendicular to the Galactic 
plane in the GC. The observed line of positions is shown in Figure 1. The emission from 
the 2-3 positions at each end of the line probably originate in the diffuse interstellar 
medium (ISM) along the line-of-sight to the GC as well as in the GC itself. Our goals are to 
clarify the source(s) of excitation of the Arched Filaments and the other thermal arcs and to 
determine the relationship of the Bubble to the clusters and other filament structures. We 
defer discussion of the continuum features seen in the spectra and the spectra of the Arches 
Cluster stars to later papers. In §2 we present the observations; in §3 we describe the results, 
including the gas excitation and atomic abundances; in §4 we present the evidence for high 
velocity shocks in the Bubble; in §5 we discuss the origins of the shocked gas in the Bubble; 
and in §6 we summarize our conclusions. 



2. Observations 

We observed the Galactic Center region with the Infrared Spectrograph (IRS) (Houck et 
al. 2004) on Spitzer Space Telescope (Werner et al. 2004) with the Short-High (SH, 10 - 19.5 
/im) and Long-High (LH, 19.5 - 38 p) modules on 21 - 23 March 2005. The LH spectra 
were taken in the IRS Stare mode, where a position on the sky is observed for a number of 
integrations ("cycles"), in this case four, and then displaced by a third the distance along 
the slit for a second group of integration cycles. Since the SH slit is half the size of the LH 
slit (4.7" x 11.3" vs. 11.1" x 22.3", respectively), the same position was mapped with the SH 
module in a 2 by 3 position map so that essentially the same area on the sky was covered. 
As with LH, four or eight integrations (three for the Arches Cluster map positions) were 
taken at each telescope pointing to provide redundancy against cosmic ray hits. The Arches 
Cluster was mapped over a ~ 45" square; however, for this paper we summed all the map 
positions so that the spectra are dominated by the interstellar gas and dust and not by the 
cluster stars (the most massive stars of the Arches Cluster are largely found in the central 
10" core, and > 80% in the central 20", Stolte et al. 2002). For each telescope pointing, the 
total integration times ranged from 18 - 120 s for SH and 24 - 56 s for LH. The positions 
and integration times were chosen such that the bright lines would not saturate the sensitive 
IRS arrays; as a result, there is no spatial agreement with positions observed with any other 
IR telescope, such as the KAO or ISO. The coordinates of each of the 38 positions are given 
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in Table 1 and are plotted in Fig. 1. 

The data were reduced and calibrated with the S13.2 pipeline at the Spitzer Science 
Center (SSC). The basic calibrated data (bed images) for each telescope pointing were 
median-combined and cleaned of rogue pixels^ and noisy order edges (the ends of the slit). 
The spectra were extracted using both the Spectroscopic Modeling, Analysis and Reduction 
Tool (SMART, Higdon et al. 2004) and the SSC extraction routines in the Spitzer IRS 
Custom Extraction tool (spiceE); we use the SMART-extracted spectra in this paper (the 
results from both spectrum extraction tools are similar). The LH spectra were defringed 
with SMART and IRSFRINGeH Because the IRS was calibrated using point sources but all 
our sources (positions) are extended, a correction for the telescope diffraction pattern as a 
function of wavelength is necessary; these "slitloss" correction factors were taken from tables 
provided with SPICE. 



2.1. Spectra and Line Fluxes 

Three examples of the observed spectra are shown in Figure 2 (corrections for offsets 
between the orders will be discussed by Simpson et al., in preparation); these positions are 
of the center of the Bubble (Fig. 2a), the W Arched Filament (Fig. 2b), and the position 
at the highest Galactic latitude (Fig. 2c). The three spectra are very representative of GC 
spectra of the Bubble, H II regions, and the diffuse ISM, respectively. 

Fluxes for the faintest (noisiest) and brightest lines were estimated by integrating over 
the line profiles (e.g., Colgan et al. 1996). These integrated fluxes are probably more 
accurate than the Gaussian fits where the line/ continuum ratio > 1 since they require no 
assumptions about the shape of the line profile; the differences with Gaussian fits are less 
than ~ 1 — 2% for these bright lines: [Ne II], [Ne III], [S III], and [Si II]. Line and continuum 
fluxes of all other lines were measured by fitting Gaussian profiles plus a linear continuum 
to each of the line profiles in the spectra using both IDL's CURVEFIT procedure and the 
procedures in SMART (Higdon et al. 2004). The Gaussian fits should be more accurate 
when the line/continuum ratio < 1, since the continuum is fit simultaneously with the line. 
Errors for all fluxes were estimated by computing the root-mean-squared deviation of the 
data from the fit; because of the systematic effects of low- level, uncorrectable rogue pixels on 



1 http:/ /ssc. spitzer. caltech.edu/irs/roguepixels 

2 http:/ /ssc. spitzer. caltech.edu/postbcd/spice. html 

3 http: / /ssc. spitzer. caltech.edu/archanaly/contributed/irsfringe 
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the spectra, these errors are consistently larger than errors estimated any other way, such as 
the statistics of median averaging or the counting statistics plus read noise from the SSC's 
pipeline. The lines discussed in this paper are listed in Table 2, along with references to 
the atomic or molecular physics parameters (wavelengths, transition probabilities, collisional 
cross sections) needed for further analysis. The measured line fluxes are given in a machine- 
readable Table 3 and selected continuum fluxes (wavelengths containing no PAH emission 
features) are given in a machine-readable Table 4, available in the on-line edition. 

2.2. Extinction 

It is well known that the GC has ~ 30 mag of visual extinction (see, e.g., Cotera et al. 
2000); this corresponds to an optical depth at 9.6 /im of rg.6 ~ 3 (Roche & Aitken 1985; 
Simpson et al. 1997). However, it is clear from looking at our extracted spectra and various 
line ratios (such as the [S III] 18.7/33.5 /im ratio) that the extinction is not constant for 
all positions. Unfortunately, because the spectra do not extend to wavelengths shorter than 
10 /im, we could not estimate r 9 . 6 by fitting extinction models to the data (e.g., Simpson 
& Rubin 1990). Instead, we estimate the extinction correction by assuming that all spectra 
have essentially the same intrinsic shape given by the ratio of the 10 and 14 /im continuum 
fluxes (where there are no strong PAH features) and fitting that ratio with a template of the 
MIR fluxes measured in the Orion Nebula (Simpson et al. 1998) in this 10 - 14 /iin spectral 
region. We will discuss this in more detail in our paper on the continuum features (Simpson 
et al., in preparation). The law for the extinction versus wavelength for the GC of Chiar 
& Tielens (2006) was used to correct line fluxes at longer wavelengths (see Appendix A). 
The extinction coefficients for each line as a fraction of the coefficient at 9.6 /im are given in 
Table 2. 



3. Results 

As can be seen in Figure 3, the strengths of the ionized lines correlate well with the 
Arches and Bubble structures seen in the radio (e.g., Fig. 1) and MIR continua (Egan et al. 
1998; Price et al. 2001). This is especially true of the lines that are strong in H II regions 
([Ne II] 12.8 /xm, [Ne III] 15.6 /im, [S III] 18.7 and 33.5 /im, the [S III] lines not being plotted 
because they are very similar to the [Ne II] line). Lines prominent in photodissociation 
regions (PDRs) ([Si II] 34.8 /im and [Fe II] 26.0 /im) also peak in the Arched Filaments 
and the Bubble rim, but there is also some contribution from the diffuse ionized gas along 
the line-of-sight to the GC. The 2-3 positions at the ends of the line almost certainly 
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sample substantial amounts of foreground gas; however, we do not subtract these line fluxes 
as "background" because (1) this diffuse ionized gas is interesting in itself and (2) the half a 
degree in Galactic latitude with the Galactic plane in the middle almost certainly precludes 
the "background" from being uniform or flat. 

3.1. Radial Velocities 

Even though the resolution of the IRS SH and LH modules is only ~ 500 km s _1 , as a 
result of our high signal/noise ratios and the IRS stability the observed radial velocities of 
the strongest lines have uncertainties as small as 10 - 20 km s _1 . These uncertainties are 
estimated from the dispersion of the differences in velocities of two lines that should have the 
same velocity at the same position, such as the two [S III] lines. These velocities are plotted 
in Figure 4. Because the wavelength calibration of the IRS is still incomplete, the velocities 
of the ionized lines were normalized to the radial velocity observed at Position 31 in the E2 
Arched Filament by Lang et al. (2001, V L sr = —15 km s _1 ). The dotted line is the average 
velocity of the four strongest ionized lines ([Si II], [Ne II], and [S III]); these velocities are 
listed in Table 1. 

The velocities of the H 2 lines were normalized as follows: The CO J = 4 - 3 intensities 
as a function of velocity and Galactic latitude were extracted from the data cube of CO 
measurements of Martin et al. (2004) at Galactic longitude 0.121°. Intensity-weighted 
average velocities were calculated for each latitude; these are plotted as a solid line in Fig. 4. 
Our H 2 velocities were normalized by requiring that the average of all our H 2 velocities be 
equal to the average of these CO velocities sampled at the same Galactic latitudes as our 
data. 

The excellent agreement of the strong, ionized lines shows that Spitzer can be used to 
measure relative radial velocities. Velocity gradients of as much as 5 - 9 km s^ 1 arcsec -1 
can even be seen in the [Ne II] and [S III] lines in a few of the SH 2 by 3 maps. The velocity 
of the weaker [Ne III] line is not as consistent — the line centroid is less well determined 
and the line may include a larger percentage of high-excitation diffuse ISM gas at a different 
velocity. 

In fact, there are indications of multiple velocity components in Bubble Positions 15, 
16, 19, 20, 21, and Bubble Rim Position 23. Not only are these the positions with velocity 
gradients in the SH ionized lines, the individual spectra for these positions have line widths 
significantly (25 - 90 km s" 1 ) larger than the instrumental line width, with Position 19 being 
the most extreme (leading to deconvolved full widths at half maximum > 250 km s" 1 ). Most 
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likely, Bubble Position 18, which has one of the highest positive velocities in Fig. 4, is missing 
the negative velocity components of the other Bubble positions and thus has relatively narrow 
lines. In general, the brighter lines are dominated by strong, single-velocity components and 
thus do not show any significant line broadening over the instrumental width. 

The observed radial velocity pattern suggests that the Arched Filaments and Bubble 
rim structure have a substantial line-of-sight component to their distances. The average 
velocity of the ionized lines (dotted in Fig. 4) shows significant structure on scales < 2 pc, 
whereas the Vlsr profile in the Arched Filaments is quite smooth. If the negative-velocity 
gas is in front of (behind) the positive-velocity gas, the negative-velocity Arched Filaments 
are closer to (farther from) the Sun than the positive-velocity, southern rim of the Bubble 
(Positions 10, 11, and 12). Position 20, with its large negative velocity, may be part of the 
northern rim of the Bubble. Some of these Bubble positions have broad lines, indicative of 
an average over both positive and negative velocity components. The fact that the Bubble 
positions (14 - 21) sometimes have broad lines and sometimes do not (but instead have large 
absolute velocities) is probably due to the ionized gas being clumpy. Positions 22 and 25 
may be the southern extension of the Sickle ("Sickle Handle" in Table 1), since they have 
positive velocity like the Sickle (e.g., Lang et al. 1997; Simpson et al. 1997). If the apparent 
Filaments and Bubble structures are caused by a wind or explosion some time in the past, 
one would expect the foreground components to have negative velocities due to expansion. 
It would be interesting to see a complete map at higher spectral resolution of the velocities 
in the Bubble region since surely, the structure seen here is affected by the sample locations. 



3.2. Molecular Hydrogen and PDRs 

Unlike the ionized lines, the H_2 line fluxes peak (Fig. 3a) near the Galactic plane 
(latitude ~ —0.05°) and have a relatively flat intensity distribution with Galactic latitude, 
similar to the distribution seen in warm CO J = 4 - 3 (Martin et al. 2004). The radial 
velocities of the H 2 lines also show a distinctly different pattern from those of the ionized 
lines (Fig. 4) (correlation coefficients < 0). This suggests that the warm molecular hydrogen 
is not associated with the same Arched Filament and Bubble structures that are seen in the 
atomic lines and dust continuum. The H_2 lines probably arise from numerous PDRs along 
the line-of-sight to the GC as well as molecular clouds in the GC itself (see the Appendix). 
We will leave additional discussion of the PDR properties of the molecular clouds to a later 
paper. 

However, we do note that the [Si II] lines (and the [Fe II] lines, although with poorer 
signal/noise) show the same velocity structure (Fig. 4) as the ionized lines (correlation co- 
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efficients > 0.8), in addition to showing a similar intensity pattern (Fig. 3) (correlation 
coefficients > 0.7). Hence, we ascribe these low-excitation ionized lines to the photoionized 
gas of the GC, rather than to the largely neutral gas associated with PDRs. 

3.3. Electron Temperatures 

Computations of H II region abundances require knowing the electron temperature, 
T e . Although the infrared forbidden lines have low sensitivity to T e (the forbidden line 
emissivities range from being proportional to T e ~ ' 25 to T~ 0A5 ), recombination lines are quite 
sensitive to T e — the emissivity of the H 7 - 6 line is proportional to T~ 13 . Measurements 
of T e have been published for the Arched Filaments by Lang et al. (2001). A typical value 
for our Arched Filament positions is T e ~ 6000 K. We adopt this value for all positions and 
include an estimated uncertainty of 500 K. In fact, we expect variations in T e and will discuss 
the effects of other values of T e in later sections. 

3.4. Electron Densities and Corrections to the Estimated Extinction 

Electron densities, N e , were estimated from the flux ratios of the density-sensitive [S III] 
18.7 and 33.5 /im lines (e.g., Rubin 1989; Simpson et al. 1995). The [S III] 18.7 to 33.5 jum 
flux ratio is also sensitive to extinction. For 8 positions the corrected 18.7/33.5 fim line ratio 
(§2.2) is below that produced by even the lowest possible N e , suggesting that an additional 
extinction correction is necessary. For these positions the estimated r 9 6 was increased by 
0.1 - 0.9 until the [S III] 18.7/33.5 jum line ratio corresponded to a minimum N e ~ 10 cm -3 . 
The estimated values of r 9 . 6 are listed in Table 1. Since both [S III] lines are very strong, 
the chief uncertainty for N e is due to uncertainties in the extinction correction. In order 
to include at least some of this extinction correction uncertainty in the error analysis, we 
assumed that the one sigma uncertainty in the optical depth is 10% of the value of rg g in 
Table 1, multiplied by the factor appropriate for the wavelength of each line (ta/t 9 . 6 ) given 
in Table 2, and propagated this uncertainty through all succeeding ionic ratio calculations. 
The computed densities are plotted in Figure 5. The highest density is at Position 31, the 
brightest observed knot in the Arched Filaments. 

Table 1 also contains the emission measure (= f N^dl) needed to produce the observed 
H 7 - 6 recombination lines. These emission measures when combined with the minimum 
N e ~ 10 cm -3 require path lengths of ionized gas ~ 100 — 500 pc. Such long path lengths 
are not unreasonable for the three end positions (1, 2, and 3) at negative Galactic latitudes 
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but they may be an indication that at some positions, especially Positions 14 and 15 in the 
Bubble, the minimum density, and hence the extinction, should be larger. For most of the 
rest of the positions, where N e > 50 cm -3 , the pathlengths are only a few pc. From this we 
infer that the gas emitting the ionized lines is clumpy, both within the Bubble and especially 
in the Arched Filaments. 



3.5. Ionic and Total Abundance Ratios 

Using the extinction-corrected fluxes and N e shown in Fig. 5, we computed various 
ionic ratios. The ratios (Ne + + Ne ++ )/H + and (S ++ + S +3 )/H + are plotted in Figure 6. 
Although we would expect the Ne/H and S/H ratios to be constant, we see that neither ratio 
is constant for all Galactic latitudes. The usual explanation for variable abundance ratios 
in a single source is that some ionization state is not being counted — this is certainly the 
case for the S/H ratio, as estimated by the (S ++ + S +3 )/H + ratio, where there could be 
significant unmeasured S + . On the other hand, for H II regions we expect > 99% of the neon 
to be either singly or doubly ionized, and even in partially-ionized diffuse ISM we expect 
neon to be essentially as ionized as hydrogen. The sizable error bars are due almost entirely 
to the measurements of the relatively weak H 7 - 6 line; the 500 K uncertainty in T e adds an 
additional ~ 10% to the error bars. We suggest that the reason for the apparent variations 
in the Ne/H abundances is our assumption of a constant T e = 6000 K. If we assumed much 
lower T e in the diffuse ISM for the positions away from the Arched Filaments, we would 
derive a much more constant (Ne + + Ne ++ )/H + ratio. For example, decreasing the assumed 
T e from 6000 K to 3000 K would increase the derived (Ne + + Ne ++ )/H + abundance ratio 
by a factor of ~ 2.1 (and the (S ++ + S +3 )/H + abundance ratio by a factor of ~ 1.9). Since 
the extinction coefficients for the H 7 - 6 and [Ne II] lines are almost identical (Table 2), the 
Ne + /H + ratio is not sensitive to extinction uncertainties. 

Because of the uncertainty in T e and the need for correction for S + in the diffuse ISM 
positions, we compute total abundances only for the Arched Filament positions. These 
positions are chosen because they should require the least correction for S + . Here the 
Si + /(Ne + + Ne ++ ) ratio is particularly low (Fig. 6c), from which we infer that silicon is 
mostly doubly ionized or more, as should be sulfur. For Positions 29 - 34 we find the aver- 
age Ne/H= 1.63 ±0.08 x 10" 4 , the average S/H= 1.16 ±0.06 x 10~ 5 , and Ne/S= 13.6 ±1.1. 
These values of Ne/H and S/H are both about 1.6 ± 0.1 times that of the Orion Nebula 
(Simpson et al. 2004, 1998). We note that with updated cross sections (Table 2; Chiar & 
Tielens 2006) and extinction (Cotera et al. 2000), the S ++ /H + ratio measured by Simpson et 
al. (1995) for G0.095±0.012, which is part of the Arched Filaments, is now 1.8 ±0.3 x 10" 5 , 
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such that S/H~ 2.0 x 10~ 5 , after correction for S+ and S+ 3 . 

Simpson & Rubin (1990), Simpson et al. (1995), Afflerbach et al. (1997), and Martin- 
Hernandez et al. (2002) showed that the gradients of abundances with respect to hydrogen 
extend all the way into the GC. Their observed gradients for Ne/H and S/H ranged from 
—0.039 to —0.086 dex kpc -1 , with a weighted average of —0.059 dex kpc -1 . However, 
the Orion Nebula has a Galactocentric distance of ~ 8.4 kpc, thus making our observed 
Arched Filament /Orion Nebula ratio of 1.6 correspond to a gradient of only —0.025 dex 
kpc -1 . Simpson et al. (1995) suggested that the gradient is flatter in the inner Galaxy 
due to the influence of the inner Galaxy bar producing radial mixing. On the other hand, 
there are variations in the chemical composition of the GC that may show that the gas is 
not well-mixed. Rodriguez-Fernandez & Martin-Pintado (2005) observed the H II region 
lines at a number of positions with the ISO spectrometers. Their measured Ne/H and S/H 
ratios (2 — 4 x 10~ 4 and 8 — 30 x 10~ 6 , respectively) are mostly larger than our measured 
ratios; however, their values include large uncertainties due to extinction and to the different 
aperture sizes at the wavelengths of Bra and the Ne and S lines. 



3.6. Sources of Excitation 

The H II region excitation can be estimated from the flux ratio of a doubly ionized 
line to a singly ionized line of the same element if the ionization potentials of both lines are 
greater than that of hydrogen (13.6 eV). For Spitzer, the best line pair is the [Ne III] 15.6 /zm 
and [Ne II] 12.8 /zm lines, since both lines have similar extinction. The excitation determined 
from the [S IV] 10.5 and [S III] 18.7 and 33.5 /im lines is less reliable because of the high 
extinction at 10 /im in the GC. The excitation ratios of Ne and S are plotted in Figure 7. 

The excitation of an H II region decreases with distance from the exciting stars, both 
because there is absorption of the high energy photons by gas closer to the stars (resulting in 
smaller Stromgren spheres for the higher ionization species) and because the radiation field 
is more dilute at larger distances. The latter effect is seen where the H II region consists of 
clumps of gas at varying distances from the star with no gas interior to absorb the high energy 
photons. This case is usually described as having a low value of the ionization parameter 
U, where U is defined as the ratio of the photon density and the electron density and is 
proportional to the photoionization rate divided by the recombination rate. The effect of a 
dilute radiation field is that the doubly and triply ionized ions tend to recombine to lower 
ionization states even if the ionizing stellar spectral energy distribution (SED) is hard enough 
that most ions are at least doubly ionized in ordinary, compact H II regions. Good examples 
are the ions of silicon and sulfur, which are normally at least doubly ionized in H II regions 
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(the singly ionized states commonly occur in PDRs). Noticeable amounts of the singly 
ionized states of these elements occur when log U < —2, regardless of the exciting SED. 
Consequently, there can be gas that is simultaneously high-excitation (containing doubly 
and triply ionized atoms) due to high-energy photons of the stellar SED (or shocks) but 
contains substantial numbers of singly-ionized atoms due to excessive recombination at low 
U. 

We can test for changes in the excitation due to distance from the exciting stars by 
observations of singly ionized silicon. The Si + /Ne ratio is shown in Fig. 6c but an even 
better test is the [Si II] 34.8/ [S III] 33.5 /tin line ratio, because it has little dependence on 
N e , T e , or extinction. Simpson et al. (1997) discussed this effect for models with varying 
U, which they compared to observations of the gas in the Sickle surrounding the Quintuplet 
Cluster. Figure 8 shows the [Si II] 34.8/[S III] 33.5 /tin line ratio for the GC positions. For 
models of compact, spherical H II regions, this ratio ranges from ~ 0.1 — 0.3 (Simpson et al. 
1997), and for observed H II regions, from ~ 0.1 to ~ 0.6 (Simpson et al. 1997; M. R. Haas, 
private communication). From Figure 8 we infer that for all the positions with the [Si II] 
34.8/[S III] 33.5 /im line ratio > 0.7, the gas lies at substantial distances from the exciting 
stars. 

We can use this apparent decrease in excitation with distance from a potential source 
of ionizing photons to help identify the actual ionizing stars. Because the Arched Filament 
Positions 29 - 34 have [Si II] 34.8/ [S III] 33.5 /im line ratios (Fig. 8) typical of compact H II 
regions, they must be ionized by stars relatively nearby. These stars must be the Arches 
Cluster stars and not the much more distant Quintuplet Cluster stars. The neon and sulfur 
excitation (Fig. 7a and 7b) decreases as one looks to higher Galactic latitude from the Arches 
Cluster as would be expected if these filaments are ionized by the Arches Cluster and are 
physically more distant. Because the gas at the position (28) of the Arches Cluster itself has 
lower excitation than the two positions at immediately higher Galactic latitudes, we believe 
that the Arches Cluster lies some distance away from the gas emitting the observed lines (see 
also Cotera et al. 2005 and Lang et al. 2001). Since the gas at Positions 22 - 26 has higher 
excitation yet (Fig. 7), we believe that this gas is excited by the Quintuplet Cluster and 
not the Arches Cluster. These positions also have radial velocities (§3.1) more in agreement 
with the Quintuplet Cluster and Sickle than with the Arched Filaments. 

We show now that the hot stars of the Quintuplet Cluster almost certainly ionize the 
rims of the Bubble. A question that motivated this project was whether the southern side of 
the Bubble rim might be ionized by hot stars lying outside the GO. 11— 0.11 molecular cloud 
to the south. We find no indication of hot stars at more southerly latitudes than the Bubble 
- such stars would produce doubly ionized Si and probably Ne. However, we observed that 
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the [Si II] 34.8/[S III] 33.5 /im line ratio is high in the rim, indicating that the exciting stars 
lie at some large distance, and the Ne ++ /Ne + ratio steadily decreases with distance from the 
Quintuplet Cluster. The region near the Sickle is probably also ionized by the Quintuplet 
Cluster with the possible exception of very low-excitation Position 20, the same position 
that has the anomalously negative radial velocity of Fig. 4. Or it may be that Position 20 is 
also ionized by the Quintuplet Cluster but has a large line-of-sight separation from it. 

We have used the photoionization code CLOUDY (version 07.02.00, last described by 
Ferland et al. 1998) to estimate the ionization structure of clumps of gas offset from the 
Quintuplet Cluster at distances ranging from 0.01 to 25 pc. We assumed empty space 
between the gas clumps and the exciting stars. The models used the 36 kK supergiant 
atmosphere of Sternberg et al. (2003), abundances equal to two times those of the Orion 
Nebula (Simpson et al. 2004), N e = 316 cm -3 , and 10 509 ionizing photons per second. The 
results for the integrated Ne ++ /(Ne + + Ne ++ ) and S +3 /(S ++ + S +3 ) ratios are plotted in 
Figure 9; also plotted are the ratios for Positions 4-26 where the abscissa is the distance of 
the position from the Quintuplet Cluster on the plane of the sky. We see that the observed 
Ne ++ /(Ne + + Ne ++ ) ratios (Fig. 9a) agree reasonably well with the models, the differences 
probably being due to the need for some additional line-of-sight contribution to the distance. 
Although the S +3 /(S ++ + S +3 ) ratios (Fig. 9b) also decrease with distance, the model ratios 
decrease much more steeply than the observations. The factor of ~ 4 discrepancy in the S +3 
fraction between the observations of positions in the Bubble rim and the models cannot be 
due to an overcorrection for extinction, since the Bubble rim has typically r 9 6 < 2, which 
produces a correction of e 2 * - 777 = 4.66 at 10.5 /im. 

The temperature of the warm dust also decreases with distance from the Quintuplet 
Cluster, as is shown in Figure 9c. Calculations of the emission from interstellar dust grains 
are usually made as a function of x, the ratio of the impinging radiation field divided by the 
"local interstellar radiation field". Li & Draine (2001) and Draine & Li (2007) calculated 
the emission spectrum of dust grains under conditions where x ranged from 0.3 to 10 4 . They 
find that the spectrum < 25 /im shows little variation with changes in the radiation field 
intensity until x > 10 3 , because for these wavelengths at the lower Xi the emission mostly 
arises from single-photon heating of very small carbon grains. For x ^ 10 3 > the 18.7/33.5 /im 
continuum ratio increases substantially with increasing x because there is now a contribution 
from larger, steady-state grains. Since x is observed to be > 10 4 in the Sickle (Simpson et 
al. 1997), it is reasonable to expect that the positions closest to the Quintuplet Cluster are 
in the regime of high enough radiation field intensity to affect the 18.7/33.5 /im ratio. The 
18.7/33.5 /xm ratio appears to flatten out at large distances from the Quintuplet Cluster 
because x is now < 10 3 - 
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Rodriguez- Fernandez et al. (2001b) showed that the excitation of both the nitrogen and 
neon lines decreases with distance from the Quintuplet Cluster by comparing their ISO line 
measurements with CLOUDY models. Most of their observations were of the Sickle region 
between the Quintuplet and Arches Clusters, but they did have five positions in the Bubble 
rim. Our observations, made with the great sensitivity of Spitzer, show that the influence of 
the Quintuplet Cluster extends well beyond the inner edge of the Bubble rim and includes 
heating the dust as well as ionizing the gas. 

4. High Velocity Shocks 

Figure 7 shows that the highest excitation occurs in the interior of the Bubble region 
and in the diffuse ISM at the ends of the line of positions. Surprisingly, the highest excitation 
does not occur close to either cluster — either there is some additional source of excitation in 
the Bubble or the gas at the positions closest to the Quintuplet Cluster actually has a large 
line-of-sight distance to the cluster. We suggest here that the source of the high excitation 
gas may be high velocity shocks in the Bubble. 

4.1. Triply Ionized Oxygen 

Because photons that can doubly ionize helium (> 54.4 eV) are required for production 
of +3 (IP = 54.9 eV), triply ionized oxygen should not be produced in Galactic H II 
regions except those ionized by the most massive O stars. However, the [O IV] 25.89 fxm line 
is detected at almost every position that we observed. The exceptions (observed signal/noise 
< 3) occur, not where the line is weak, but where the continuum, with its noise-producing 
rogue pixels, is very strong. The normalized line profiles that show this line are plotted in 
Figure 10. Fluxes in the LH beam (248 square arcsec) with > 3<r measurements range from 
2.1 x 10~ 17 W m -2 to 1.5 x 10~ 16 W m -2 . The peak fluxes occur in Positions 19 - 27 but 
the best signal/noise is for positions in the Bubble and the diffuse ISM at high Galactic 
latitudes where the line/continuum ratio is the largest. We plot the +3 /(Ne + + Ne ++ ) 
ratio in Figure 11. Since the O/Ne ratio, which is expected to be approximately constant, 
is observed to equal ~ 4 in Galactic H II regions and +3 /(Ne + + Ne ++ ) is < 0.01 here, we 
see that very little of the oxygen is +3 and therefore very little neon should be in higher 
ionization stages, too. 

Although triply ionized oxygen should not be produced in Galactic H II regions, it has, in 
fact, been detected in some extra-galactic H II regions. However, either these are extremely 
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low metallicity, which hardens the stellar ionizing SED, or the exciting stars include hot 
WR stars (Schaerer & Stasihska 1999). This is not the case in the GC — we observe no 
enhancement of the +3 / (Ne + + Ne ++ ) ratio at the position of the Arches Cluster, indicating 
that its WR stars are not effective at producing +3 . 

Lutz et al. (1998) discussed the [O IV] that they had observed with ISO in some galactic 
nuclei that are not low metallicity. They attempted to reproduce the [Ne III]/[Ne II] and 
[O IV]/([Ne II] + 0.44[Ne III]) line ratios using CLOUDY (Ferland et al. 1998). They 
modeled both a weak active galactic nucleus embedded in a starburst and an H II region 
ionized by main-sequence stars plus a hot (80 kK) blackbody to represent WR stars. All 
their models showed the same pattern: in order to produce as much [O IV] emission as 
is observed, they had to include so much of the hot component that the [Ne III]/[Ne II] 
ratio was much larger than observed. They noted that if there were small regions producing 
relatively strong [O IV] intermixed with low excitation gas, they might be able to produce the 
observed ratios; however, they considered it unlikely because they knew of no such regions 
in the Galaxy or 30 Dor, and planetary nebulae produce too small a contribution to be 
detected. They suggested that the high excitation of [O IV] is caused by shocks. 

We have used most recent version of CLOUDY (07.02.00) to make similar models, 
employing both black bodies and the supergiant SEDs of Sternberg et al. (2003), who used 
the WM-BASIC model atmosphere code of Pauldrach et al. (2001). Supergiant models are 
defined here as the atmosphere models with the smallest gravity, g, for each stellar effective 
temperature, T e g. Supergiant SEDs are used because H II region models computed with these 
SEDs as input give substantially better agreement with observed Ne ++ lines than the WM- 
BASIC SEDs with other gravities or computed with other model atmosphere codes (Simpson 
et al. 2004; Rubin et al. 2007). We emphasize that even though the actual stars may be 
something other than supergiants, the integrated cluster SEDs are similar to the calculated 
supergiant model SEDs. The CLOUDY models have as input density = 316 cm~ 3 ; Sternberg 
et al. SEDs with T e s = 32 kK to 50 kK; number of ionizing photons A^L yc = 10 49 , 10 50 , and 
10 51 s -1 ; and inner shell radius Ri= 0.01, 1.0, and 10 pc. The blackbody models have T c g = 
30 kK to 45 kK and similar other parameters. Many of the CLOUDY models with blackbody 
SEDs have additional blackbody components with T ef[ = 10 5 and/or 10 6 K to represent hot 
WR stars or the pervasive X-ray flux of the GC. 

Figure 12 shows the fraction of +3 versus the fraction of Ne ++ for both the H II region 
models and our observations. The only models using Sternberg et al. (2003) atmospheres 
(open circles) with enough +3 to show up on this plot are the two models with T cS = 50 
kK, N Lyc = 10 49 s" 1 , and Ri = 0.01 or 1.0 pc (since star clusters containing multiple 50 kK 
stars are not likely in the Galaxy, no models with 50 kK atmospheres and larger A"L yc were 
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computed). The asterisks are the blackbodies with additional hot components for WR stars 
and X-rays. As was earlier shown by Lutz et al. (1998), all the data fall to the upper left 
(larger +3 /Ne ++ ) with respect to all the models calculated with believable SEDs (asterisks 
and open circles). 

Before leaving this topic we test whether there is any SED that could be used in an 
H II region model that would reproduce the observed ionic ratios. To do this we computed 
models with the Sternberg et al. (2003) SED for 36 kK, (which has a sharp drop-off at 54.4 
eV), but we arbitrarily increased the flux for photon energies E > 54.5 eV (> 4.013 ryd). 
The other inputs were Ri = 10 pc and log A^yc = 50.5. Such SEDs do not affect ionic ratios 
of low-excitation ions, unlike SEDs where the flux at 54.5 eV is increased by the addition of a 
blackbody, which has substantial photon emission at energies lower than its maximum. The 
ionic ratios from these models are plotted as open triangles in Fig. 12. The lowest triangle 
on this plot is from the model that has a jump in the SED of ~ 4 orders of magnitude going 
from 54 eV to 54.5 eV. In this model the SED is increased by a factor of 10 9 over the original 
very low flux with E > 54.5 eV found in Sternberg et al's model atmospheres. The open 
triangles connected by the solid line have an increase of a factor of 10 10 over the original 
SED; their values of Ri are 2, 5, 10, and 15 pc (top to bottom). We see that the only model 
lying in the area of the observations is the model with Ri = 10 pc and a SED with a flux 
jump of a factor of ~ 10 5 at 54.5 eV. SEDs with other flux jumps and Ri could probably 
fit, but if the SED represents the ionizing source of all the positions, which all have different 
Ri, it would not be possible for the same SED to fit all positions. Moreover, it is difficult to 
imagine a physical process that would produce such strong emission in the He ++ continuum 
> 54.5 eV when the normal stellar SED has deep absorption. We conclude that the +3 is 
not produced by photoionization. 

The advantage of observing [O IV] in the GC is that we can locate the maximum of the 
high-excitation gas as being in a region that does not appear to contain the exciting stars. 
We also can eliminate any possibility that the [O IV] is confined to small but high excitation 
regions as we detected it at all our positions. According to the shock models of Contini & 
Viegas (2001), a shock velocity ~ 100 km s" 1 is adequate to produce the amount of [O IV] 
that is observed. 



4.2. Gas-Phase Iron Production in Shocks 

Whereas +3 could be produced by high energy photons such as X-rays, high gas-phase 
abundances of elements normally depleted onto grains indicate grain destruction in shocks. 
Here we demonstrate that the gas-phase iron abundance is greatly increased in the center of 



-17- 



the Bubble compared to the surrounding gas. This is suggestive of shock production. 

To compute the Fe abundance ratios we modeled the Fe + and Fe ++ ions with 36 and 
34 level statistical equilibrium calculations, respectively (see Table 2 for atomic parameters 
and references). Figure 10 shows the observed [Fe II] lines and Figure 13 shows the derived 
(Fe + + Fe ++ )/(Ne + + Ne ++ ) ratio. Here we have plotted ratios computed using both the 
[Fe II] 25.99 jum line and the [Fe II] 17.94 /iin line since the former includes substantial 
emission from the low temperature PDRs along the line-of-sight (e.g., Kaufman et al. 2006) 
but the lower level of the latter is 1873 cm" 1 above ground and the line should arise only in 
the H II region gas, like the [Fe III], [Ne II], and [Ne III] lines. Unfortunately, the [Fe II] 17.94 
jum line is weak and is somewhat blended with the [P III] 17.885 /iin line (Fig. 10a). With 
the Spitzer resolution of ~ 600, the lines are separated by ~ 1.7 resolution elements and so 
can usually be distinguished; however, the measurements have some additional uncertainty 
because the 2-line- Gaussian-fit routines in SMART would converge only if the line width was 
fixed at the Spitzer IRS SH resolution. There are two other [Fe II] lines in the LH bandpass 
but they are either in a very noisy part of the wavelength range (35.35 /im) or in a region 
with bad fringing (24.52 /mi). The latter line is usually detectable, but since its signal/noise 
ratio is much worse than the 17.94 /im line, we do not use it. The [Fe III] 33.04 /tin line also 
is usually detectable (Fig. 2a), but since it imparts no additional information with respect to 
N e or T e and its signal/noise ratio is so much poorer than the 22.9 /iin line that the 33.0/22.9 
line ratio is not usable for extinction estimates, we do not use it either. 

We estimate the gas-phase abundances of iron in the GC by using our observed (Fe + + 
Fe ++ )/(Ne+ + Ne ++ ) ratios (Fig. 13): (Fe + + Fe ++ )/(Ne+ + Ne ++ )= 0.0066 ±0.0016 in the 
Arched Filaments and (Fe+ + Fe ++ )/(Ne+ + Ne++)= 0.043 ± 0.004 in the Bubble. Since 
Fe ++ is moderately-low excitation, there can be substantial Fe +3 in photoionized H II regions 
(see Rodriguez & Rubin 2005 for a discussion of the uncertainties in the measurements of Fe 
abundances in H II regions). We estimate the correction for Fe +3 using both the CLOUDY 
models described above and models computed with the code NEBULA (Simpson et al. 
2004; Rubin et al. 2007), noting that NEBULA and CLOUDY produce similar results for 
the same SEDs and other input parameters. The ionization fractions computed with these 
models are plotted in Figure 14. As explained above, we prefer the models computed with 
the supergiant atmospheres of Sternberg et al. (2003), which are connected with the heavy 
solid line in the plot; from these we estimate an ionization correction factor of ~ 20%. This 
gives the gas-phase (Fe/Ne) gas = 0.0082 and 0.054 in the Arched Filaments and Bubble, 
respectively. Assuming that the Ne/H ratio observed in the Arched Filaments (Ne/H = 
1.63 x 10~ 4 ) is applicable to all the GC gas and that Ne is not condensed onto grains, by 
multiplying (Fe/Ne) gas by Ne/H we estimate that (Fe/H) gas ~ 1.3 x 10 -6 in the Arched 
Filaments and ~ 8.8 x 10~ 6 in the Bubble. There is an increase in the gas-phase abundance 
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of Fe of a factor of 6.5 ± 1.7 in the Bubble with respect to the surrounding gas. 

In order to estimate the depletion factors for Fe in the GC, one first needs to know the 
total Fe abundance in both grains and gas. There are several ways to do this: 

(1) There are interstellar absorption lines of O, Ne, and Fe from 13.4 - 23.5 A that can 
be observed by Chandra and XMM-Newton. The advantage of X-ray measurements is that 
the atoms in the grains produce sufficient absorption that total abundances are measured, 
not just gas phase abundances. Using Chandra, Juett et al. (2006) observed absorption lines 
and edges in front of 7 neutron stars, measuring averages of O/Ne = 5.4 ± 1.6 and Fe/Ne 
= 0.20 ± 0.03 (unfortunately, this soft X-ray wavelength range suffers too much interstellar 
extinction for observations of GC sources). They derive the expression 

(Fe/Ne)isM = (Fe/Ne) meas /(1 -0 + fP) (1) 

to correct for optical depth effects in the larger grains, where (Fe/Ne)isM is the total gas 
plus grains ratio, (Fe/Ne) meas is the measured average ratio given above, f3 is the fraction of 
Fe residing in grains, and / is the factor accounting for grain optical depths and the grain 
size distribution. They derive / = 0.88 for Fe and / = 0.90 for O. Using eq. (1) we find that 
(Fe/Ne) ism ranges from 0.20 to 0.23 and (0/Ne)i SM ranges from 5.4 to 6.1 for < p < 1. 

So that they would not need to make corrections for grains, Yao et al. (2006) measured 
high ionization absorption lines of Ne and Fe in the low-mass X-ray binary 4U 1820—303, 
with the results that O/Ne = 6.2 and Fe/Ne = 0.30. 

(2) There are X-ray emission lines of Fe at 6.4 and 6.7 keV. Wang et al. (2006) measured 
the spectra of the Arches Cluster stars with Chandra. From the 6.7 keV line — which probes 
highly ionized gas in the shocked cluster winds — they estimated that the Fe abundance 
is 1.8 times Solar, referring to the Solar abundances of Anders & Grevesse (1989), or Fe/H 
= 8.4 x 10" 5 . From measurements of the spectrum of the diffuse emission surrounding the 
Arches Cluster with Suzaku, Tsujimoto et al. (2007) estimated that the abundance of Fe 
is Solar, again with respect to the abundance summary of Anders & Grevesse, such that 
Fe/H~ 5 x 10~ 5 . Koyama et al. (2007), using Suzaku, estimated an average Fe column 
density= 9.7 x 10 18 cm" 2 for the diffuse gas in the GC, from which they calculated that 
the iron abundance is Fe/H~ 1.6 x 10~ 4 , assuming that the column density of H is 6 x 10 22 
cm" 2 . For these iron abundances, we estimate Fe/Ne~ 0.5, 0.3, or 1.0, respectively. 

(3) By observing the spectra of cool stars in the GC, Ramirez et al. (2000) and Blum 
et al. (2003) found that the abundance of iron in the GC central cluster is ~ Solar. From 
the spectra of hot stars in the central cluster and the Quintuplet cluster, Najarro (2006) 
also found a Solar abundance of iron. These results are in contradiction to the computations 
giving a steep gradient for Fe in the Galaxy (e.g., Cescutti et al. 2007). The most recent 
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compilation of Solar abundances is that of Asplund et al. (2005), who give Fe/H= 2.8 x 10 5 . 
From this we estimate that the total Fe/Ne~ 0.17. 

In summary, estimates of the total Fe/Ne ratio in the GC range from ~ 0.17 to ~ 1.0. 
With these estimates (0.17/0.30/1.0) of the Fe/Ne ratio and our measurements, we find that 
the fraction of Fe in the gas phase is 0.048/0.027/0.0082 and 0.32/0.18/0.054 for the Arched 
Filaments and the Bubble, respectively. Thus, in spite of the factor of 6.5 increase of gas- 
phase iron in the Bubble with respect to the surrounding ISM, most Fe is still incorporated 
in grains and is not in the gas-phase. 

Dust grains can be destroyed in shocks, both by grains shattering into smaller particles 
and being vaporized in grain-grain collisions and by removal of individual atoms by photons 
and by high- velocity gas (see Jones 2004 for a review). In high velocity shocks, V s > 200 
km s _1 , the grains are stopped in the hot gas. Since in hot gas H + and He + have the 
same kinetic energies (~ \kT) but H + is more abundant, sputtering by protons dominates. 
In slower shocks, however, most of the destruction is done by the betatron acceleration of 
grains, which then move at high speeds (~ V s ) through the gas. For this case, the collision 
velocity is the same for H + and He ++ , and since the He ions have 4 times the energy of the 
protons, they dominate the sputtering, even though they are less abundant by a factor of 
ten. 

Jones et al. (1994, 1996) computed grain destruction models for a variety of shock 
velocities, initial densities n , and magnetic fields B . Both n and B in the GC are 
probably higher than those used in any of their models (especially B ) but the destruction 

— 1/2 

rate should scale as Bon . The observed depletions in the Bubble correspond to shock 
velocities ~ 100 — 50 km s -1 , depending on which Fe gas-phase fraction is correct. These 
velocities are plausible if the shocks are due to the expansion of the Bubble (see the range 
of velocities in Table 1). 

We note that there is no obvious increase in the Si + /H + ratio in the Bubble compared to 
the Bubble exterior that cannot be attributed to normal H II region photoionization, where 
most Si is Si ++ or higher. In the Bubble the average Si + /(Ne + +Ne ++ )= 0.102 ± 0.003, from 
which we estimate that Si + /H~ 1.7 x 10~ 5 . This is about 0.5 times Solar (e.g., Asplund et 
al. 2005). However, if one of the higher depletion estimates for Fe is correct, it would require 
only a small fraction of the grains to be destroyed in the Bubble, producing a non-detectable 
increase in the gas-phase Si abundance. 
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5. Discussion 

Bubbles around hot stars and star clusters are commonly observed in other galaxies 
(e.g., Oey & Massey 1994; Oey 1996) and in the Milky Way (Churchwell et al. 2006). Early 
analytic models were made by Weaver et al. (1977), who modeled the time- varying stellar 
wind of a main-sequence star expanding into the ISM. Later, Garcia-Segura & Mac Low 
(1995a, 1995b) made both analytic and numerical models of stellar winds from WR stars. 
When these winds encounter gas from previous mass-loss episodes (main- sequence and red 
super-giant stages) and/or interstellar gas, a shock results, producing high temperatures and 
hence high pressures. A shell is formed surrounding the high-pressure interior; the expanding 
shell sweeps up the exterior interstellar gas, shocking and heating it. The temperatures 
of the lower-density bubble interior are sufficiently high that, at least sometimes, X-rays 
are produced (Chu & Mac Low 1990). Although these models do predict the formation of 
bubbles, Oey (1996) found that the predicted bubbles are much larger than observed. Oey & 
Garcia-Segura (2004) found that numerical hydrodynamic and radiative models that include 
substantial pressure, P/k ~ 10 5 cm -3 K, in the ambient ISM result in bubbles and shells 
much more similar in size to the observed bubbles. 

The Bubble seen in the GC is fairly small (radius ~ 10 pc) compared to the bubbles 
in the Large Magellanic Cloud observed by Oey (1996) and Oey & Garcia-Segura (2004). 
On the other hand, the ambient ISM pressure in the GC is probably at least an order of 
magnitude larger — for the GC as a whole, Spergel & Blitz (1992) estimate P/k ~ 5 x 10 6 
cm -3 K. Rodriguez- Fernandez et al. (2001a) measured CO lines with the IRAM 30 m 
telescope and H 2 lines with ISO at their position M+0.16— 0.10 in the center of the Bubble. 
They estimated an H 2 density of ~ 10 4 cm -3 from their CO measurements and a temperature 
of 157 — 900 K from their H 2 line measurements, giving P/k = nT ~ few xlO 6 cm -3 K. 
This gas is most likely found in a PDR, probably the edge of the molecular clouds in this 
region. The turbulent pressure must be large, given that the H 2 line widths (180 km s _1 ) 
are much larger than the sound speed. Possibly the largest contributor to the ISM pressure 
is the magnetic field, B. Yusef-Zadeh & Morris (1987a) suggested that the magnetic field 
in the linear non-thermal filaments of the Radio Arc is as high as 1 mG. If so, the magnetic 
pressure (Pb/A; = 2.9 x 10 8 -B^ tG cm -3 K) dominates the gas pressure, although it would be 
found in both the interior and the exterior of the Bubble. If the linear filaments do not 
thread through the Bubble, B could be much lower. LaRosa et al. (2005) found that the 
global field of the GC is only of order 10 fj,G with an upper limit of < 100 /iG. Since the 
magnetic pressure is proportional to B 2 , the global magnetic pressure, in this case, would 
not contribute significantly to the ambient ISM pressure. 

Looking at the models of Oey & Garcia-Segura (2004), we suspect that similar models 
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with ambient ISM pressure > 10 6 cm -3 K could produce bubbles with radii of order 10 - 
20 pc, such as are observed in the high-pressure environment of the GC. The star clusters 
producing the bubbles in their models have ionizing luminosities A^yc = 10 49 ' 4 — 10 50 2 s~ x 
(we include here only their "Pre-SN" clusters because these bubbles have thicker observed 
shells like the Bubble and not the very thin shells of "Post-SN" clusters). The Quintuplet 
Cluster lies within the Bubble, although not at its center, and has iV Lyc ~ io 50 - 9 photons s _1 
and M > 7 x 10 -4 M yr _1 (Lang et al. 2005). Since the sizes of observed expanding shells 
are always much smaller than what is expected from the kinetic energy of the stellar wind, we 
speculate that the Bubble is produced by the Quintuplet Cluster. The Bubble asymmetry 
is probably caused by the asymmetry in the surrounding molecular clouds, such that the 
Quintuplet Cluster is close to the cloud next to the Sickle (Simpson et al. 1997). Before the 
Cluster wind started creating the Bubble, there was a gradient in the gas extending from 
the Vlsr = +25 km s -1 GO. 20— 0.033 molecular cloud (Serabyn & Glisten 1991) through the 
location of the stars to lower Galactic latitudes. 

We can ask about the extent of the shocked gas in the Galactic longitude direction. 
The Spitzer archive includes Program 18 (PI: J. Houck), which contains SH and LH spectra 
(all short integrations) of 18 FIR sources from Galactic longitudes —0.96 to +3.06. We 
searched these data for the [O IV] line and other indicators of high-velocity shocked gas. 
There is one position in the middle of the Bubble with [O IV] at similar intensity to our own 
Bubble positions; the two positions on either side, at Galactic longitude —0.15 and +0.24, 
have [O IV] line fluxes factors of 4 - 6 fainter. The next positions on either side have no 
[O IV]. The [Fe III] flux also peaks in their Bubble position. The relatively small extent of 
the shocked gas in Galactic longitude is another indication that the Bubble is local to the 
Quintuplet Cluster region. 

Rodriguez- Fernandez & Mart in- Pint ado (2005) observed most of these positions with 
ISO. They also found the gas in the center of the Bubble (position M+0.16— 0.10) has the 
highest excitation as seen by the [N III]57/[N II]122 /im and [Ne III] 15.6/[Ne II] 12.8 /iin 
line ratios. They conclude that the H II region gas at this position is mainly excited by 
the Quintuplet Cluster but the gas seen in lines-of-sight more distant from the Arches and 
Quintuplet Clusters must have their own ionizing sources since their ionization does not 
decrease with distance from the clusters. This is in accord with our conclusion that the 
Bubble is excited by the Quintuplet Cluster stars. 
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6. Summary and Conclusions 

We present Spitzer IRS spectra of the Arched Filaments, Sickle Handle, and Bubble in 
the GC, along with spectra of the diffuse ISM at Galactic latitudes up to ±0.3° from the 
Galactic plane. In this paper we concentrate on interpreting the ionic lines, using the H 2 
lines and continuum mainly for estimation of the interstellar extinction. 

From our estimates of the radial velocities of the H2 lines, accurate to 10 - 20 km s _1 , 
we infer that there are significant contributions to the H2 lines seen in our spectra from 
gas at locations along the line-of-sight other than the location of the ionized gas. The H 2 
line ratios are consistent with emission from multiple components at different temperatures. 
The low excitation ionized lines ([Si II] and [Fe II]) have both radial velocities and relative 
intensities better correlated with those of the more highly ionized H II region lines than with 
the H2 lines. We conclude that these lines are formed in both the H II regions and in PDRs 
closely associated with these H II regions. 

The abundances of neon and sulfur with respect to hydrogen are a factor of ~ 1.6 larger 
than those of the Orion Nebula, the standard for interstellar abundance studies. We note 
that these abundances are quite sensitive to T e because of the sensitivity of the H 7 - 6 line 
emissivity to T e , which we assume equals 6000 ±500 K. The effect is such that if one assumes 
a lower T e , as though there were more cooling from higher metallicity, one will as a result 
deduce higher metallicities. This shows the importance of accurate measurements of T e . 

The measurement positions that we chose sample both sides of various ionization struc- 
tures, thereby enabling identification of the direction from which the ionizing photons origi- 
nate. The south rim (b ~ —0.17° to —0.25°) of the Bubble is ionized by a source at higher 
Galactic latitude, probably the Quintuplet Cluster. This part of the Bubble rim is not the 
outer edge of the GO. 11— 0.11 molecular cloud, which lies along the same line-of-sight. 

The Arched Filaments are ionized by the Arches Cluster. The Arched Filament positions 
are the only locations where the Si + /H + abundance ratio is relatively low, probably because 
silicon is largely doubly ionized there. Such an ionization pattern can happen only when the 
ionizing stars are nearby, and the Arches Cluster is much closer to the Arched Filaments 
than any other known source that produces a substantial number of ionizing photons. 

The highest excitation, as seen by the indicators Ne ++ /Ne + , S +3 /S ++ , Si + /S ++ , and 
+3 /(Ne + + Ne ++ ), is found in the center of the Bubble and not at the locations of the 
Arches and Quintuplet Clusters. At the present time there is probably a several pc line-of- 
sight offset between the Arches Cluster and the ionized gas. This is not unreasonable since 
the clusters are of ages such that they have orbited the Galaxy at least once since their 
formation. 
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There are multiple indicators of strong shocks, peaking in the center of the Bubble. 
This is indicated both by the presence of +3 , which cannot be photoionized by stars with 
as large a metallicity as exists in the GC, and by the increase by a factor of 5 - 8 of gas-phase 
iron from Fe/H ~ 1.3 x 10~ 6 in the Arched Filaments to Fe/H ~ 9 x 10~ 6 in the Bubble. 
This increase in the gas-phase iron abundance is probably produced by shock destruction of 
grains. The strong shocks possibly indicate the location of a violent outflow in the recent 
past. We suggest that the source of this outflow is the Quintuplet Cluster. A shock velocity 
< 100 km s" 1 is adequate to produce both the observed +3 and Fe gas . 

Finally, we note that we have observed a wide variety of conditions in the GC. These 
would all occur in the same aperture in an observation of a distant starburst galaxy. One 
cannot assume that the highest excitation gas is associated with the most massive cluster 
— in the GC it most certainly is not. 

This work is based on observations made with the Spitzer Space Telescope, which is op- 
erated by the Jet Propulsion Laboratory, California Institute of Technology, under a contract 
with NASA. Support for this work was provided by NASA. JPS acknowledges support from 
NASA/Ames Research Center Research Interchange Grant NNA05CS33A to the SETI In- 
stitute. RHR acknowledges NASA support through LTSA and Spitzer grants. The IRS was 
a collaborative venture between Cornell University and Ball Aerospace Corporation funded 
by NASA through the Jet Propulsion Laboratory and Ames Research Center. SMART was 
developed by the IRS Team at Cornell University and is available through the Spitzer Science 
Center at Caltech. We thank M. R. Haas and A. G. G. M. Tielens for reading the manuscript 
and we thank the referee for the thoughtful and thorough review which improved the paper 
significantly. 

Facility: Spitzer (IRS). 

A. Extinction Laws 

The extinction law of Chiar & Tielens (2006) for the GC is notable for having a much 
larger 18 /im silicate feature (with respect to the 10 /im silicate feature) than the older 
relations of Draine (Draine & Lee 1984; Draine 2003a, b). With the extinction as high as it 
is in the GC, the choice of extinction law for the extinction correction obviously affects the 
various MIR line ratios. 

We can test this effect using the three observed H 2 lines. Figure 15 shows flux ratios of 
the extinction-corrected H 2 S(2)/S(l) and S(0)/S(2) lines plotted versus the optical depth at 
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9.6 /iin, r 9 . 6 (the 10 jum silicate feature peaks at 9.6 fim in Chiar & Tielens' (2006) table for 
the GC). Because Draine's (2003a, b) relation is substantially different, we also estimate r 9 6 
for each position using his law for the diffuse ISM (R = 3.1) and then correct our observed 
H 2 line fluxes for this extinction. The ratios of the fluxes thus corrected are also plotted in 
Fig. 15. We see that the H2 flux ratios corrected with Draine's extinction law have a strong 
dependence on r 96 whereas the flux ratios corrected with Chiar & Tielens' GC law have a 
much smaller dependence on r 9 . 6 . Since the extinction at 28 /im is due to the tail of the 18 
jum silicate feature, the S(0)/S(1) flux ratio does not show any difference due to extinction 
law and is not plotted; the slope is slightly negative, perhaps indicating that the positions 
with higher extinction (mostly in or near the Arched Filaments) have warmer H 2 . 

Figures 16 and 17 shows the H 2 S(0)/S(1) and S(0)/S(2) flux ratios plotted versus the 
S(2)/S(l) flux ratio for our observed fluxes corrected for extinction with the laws from both 
Chiar & Tielens (2006) for the GC and from Draine (2003a, b). We also plot curves of the 
theoretical equilibrium ratios for local thermodynamic equilibrium (LTE). The observed flux 
ratios fall far above this line; in effect, the S(l) line is much weaker than would be predicted 
from the S(0)/S(2) ratio. The LTE curve assumes the equilibrium ortho/para ratio (an 
extra factor of 3 in the statistical weights for the odd levels) — deviations such as seen here 
have been interpreted as non-equilibrium excitation due to shocks. However, it is not likely 
that the observed H 2 emission comes from a region of constant physical conditions along the 
line-of-sight to the GC. The dashed and dotted lines show that the observed ratios can be 
produced by simple two-component models: a cooler molecular cloud that produces most of 
the H 2 S(0) flux and a warmer PDR that produces most of the H 2 S(2) flux. However, we 
see that the flux ratios corrected with the Draine (2003a, b) extinction law require a much 
larger high-temperature component than the flux ratios corrected with the Chiar & Tielens 
(2006) extinction law. We conclude from Figs. 15 to 17 that the Chiar & Tielens (2006) GC 
extinction law gives a better fit to the H 2 line flux ratios than does the Draine (2003a, b) 
extinction law. 
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Table 1. Observed Positions 
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29 


17 45 47.4 


-28 48 35 


0.1274 


0.0340 


3.11 


-13.2 


2.3E5 




30 


17 45 46.5 


-28 47 59 


0.1342 


0.0420 


3.12 


-19.3 


3.3E5 




31 


17 45 45.0 


-28 47 49 


0.1337 


0.0481 


3.40 


-15.0 


7.2E5 


E2 Filament 


32 


17 45 40.7 


-28 48 04 


0.1219 


0.0595 


2.51 


-19.9 


1.7E5 




33 


17 45 35.5 


-28 47 25 


0.1214 


0.0812 


3.21 


-18.3 


2.5E5 


W Filament 


34 


17 45 34.0 


-28 47 15 


0.1209 


0.0873 


3.02 


-21.7 


5.0E5 


W Filament 


35 


17 45 31.9 


-28 46 50 


0.1228 


0.0975 


2.00 


-13.8 


2.3E5 


W Filament 


36 


17 45 21.8 


-28 45 26 


0.1235 


0.1411 


0.80 


-6.9 


1.1E4 


Diffuse ISM 


37 


17 45 10.5 


-28 44 04 


0.1214 


0.1882 


0.42 


-7.5 


1.2E4 


Diffuse ISM 


38 


17 44 58.8 


-28 42 56 


0.1152 


0.2345 


0.33 


-7.0 


9.6E3 


Diffuse ISM 



a Positions whose values of T9.6 were increased to produce the minimum [S III] 18.7/33.5 (im line ratio after correction 
for extinction. 

b This maser source (Glisten & Downes 1983) has the radio appearance of an Ultra-Compact H II region. Our position 
is slightly offset from the radio source to avoid foreground stars. 
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Table 2. Line Parameters 



Line 




Wavelength 
(fim) 


jpa 

(eV) 


Pupper 

(cm- 1 ) 


T\/T 9 . 6 


References 


TT o / r\\ 

H 2 S(0) 




28.219 





O f A O T 

354.37 


0.369 


17, 18 


t t n/i \ 

H 2 S(l) 




17.035 





705.69 


0.502 


17, 18 


TT O /o\ 

H 2 S(2) 




12.279 





■1 -1 n n 

1168.80 


1 \ f\ A A 

0.344 


17, 18 


TT T '"7 Z" 1 

H I 7-6 




12.372 


13.60 


107440.45 


0.332 


1 

13 


I / \ X T 7-"] 2t~» 

[0 IV] z P 3/2 - 


2n 
-Pl/2 


or or\r\ 

25.890 


54.94 


386.245 


1 \ A 1 \ A 

0.404 


A 1 r 

4,15 


Tat ttI 2ti 

[Ne II] 2 Pi/ 2 


- P3/2 


12.814 


"i r r 1 

21.56 


780.42 


C\ O "1 

0.313 


12 


Tat tttI .' 1 t ^ 

[Ne IIIJ Fi- 


" P2 


1 V P" P" P" 

15.555 


40.96 


642.88 


0.383 


a n 

4, 6 


re* ttI 2t~» 

[Si II] 2 P 3 / 2 - 


2n 
"Pl/2 


/i 01 r 

34.815 


1 r 

8. 15 


287.23 


0.306 


O A 

2, 4 


[P IIIJ ^P 3 /2" 


2n 
-Pl/2 


17.885 


19.77 


559.13 


0.539 


11 


[S III] 3 Pi- 3 


Pn 

U 


33.481 


23.34 


298.68 


0.311 


4, 16 


[S III] 3 P 2 - 3 


Pi 


18.713 


23.34 


833.06 


0.548 


5, 16 


[S IV] 2 P 3/2 - 


" 2 Pl/2 


10.511 


34.79 


951.43 


0.777 


11 


[Cl II] 3 Px- ? 


P 2 


14.369 


12.97 


696.00 


0.310 


14 


[Fe II] 6 D 7/2 - 


- 6 D 9/2 


25.988 


7.90 


384.79 


0.403 


1, 7, 9, 10 


[Fe II] 6 D 5/2 - 


- 6 D 7/2 


35.349 


7.90 


667.68 


0.305 


1, 7, 9, 10 


[Fe II] 4 F 7/2 - 


- 4 Fg/2 


17.936 


7.90 


2430.10 


0.540 


1, 7, 9, 10 


[Fe II] 4 F 5/2 - 


- 4 F7/2 


24.519 


7.90 


2837.95 


0.425 


1, 7, 9, 10 


[Fe III] 5 D 3 - 


- 5 D 4 


22.925 


16.19 


436.21 


0.448 


3, 8, 19 


[Fe III] 5 D 2 - 


5 D 3 


33.038 


16.19 


738.88 


0.314 


4, 8, 19 



ionization Potential (IP) is the energy required to produce the ground 
state of the molecule or ion. 



References. — (1) Bautista & Pradhan 1996; (2) Dufton & Kingston 1991; 
(3) Erickson et al. 1989; (4) Feuchtgruber et al. 1997; (5) Kelly k Lacy 1995; 
(6) McLaughlin & Bell 2000; (7) Nussbaumer & Storey 1988; (8) Quinet 1996; 
(9) Quinet et al. 1996; (10) Ramsbottom et al. 2005; (11) Saraph & Storey 
1999; (12) Saraph & Tully 1994; (13) Storey & Hummer 1995 (14) Tayal 2004; 
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(15) Tayal 2006; (16) Tayal & Gupta 1999; (17) Turner et al. 1977; (18) 
Wolniewicz et al. 1998; (19) Zhang 1996. 
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Table 3. Observed Line Fluxes a 



Position 


Galactic 
Longitude 


Galactic 
Latitude 


Flux(H 2 S(0)) b Error(H 2 S(0)) b Flux(H 2 S(l)), etc. b 
(W m" 2 ) (W m" 2 ) (W m" 2 ) 


1 
2 


0.1233 
0.1028 


-0.2998 
-0.2837 


2.076e-16 3.707e-18 5.787e-16 
2.246e-16 2.661e-18 4.854e-16 



a Table 3 is published in its entirety in the electronic edition of the Astrophysical Journal. 
A portion is shown here for guidance regarding its form and content. 

b The electronic edition of this table contains the observed fluxes and errors per LH 
beam area of the H 2 S(0), S(l), and S(2) 28.92, 17.04, and 12.28 fim lines, the H 7 - 6 
12.37 /im line, the [O IV] 25.89 /im line, the [Ne II] 12.81 /im line, the [Ne III] 15.56 /im 
line, the [Si II] 34.82 /im line, the [P III] 17.89 /xm line, the [S III] 33.48 and 18.71 /im 
lines, the [S IV] 10.51 /im line, the [Fe II] 25.99 and 17.94 /im lines, and the [Fe III] 22.93 
/im line. 



Table 4. Observed Continuum Fluxes a 



Position 


Galactic 


Galactic 


Flux(10.0-10.48 ^m) b 


Flux(13. 5-14.3 /im) b 


Flux(15.56 ^m) b 


Error(15.56 fim), ctc. b 




Longitude 


Latitude 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


1 


0.1233 


-0.2998 


2.32E-01 


6.39E-01 


5.36E-01 


5.60E-03 


2 


0.1028 


-0.2837 


2.11E-01 


5.86E-01 


4.81E-01 


4.85E-03 



a Table 4 is published in its entirety in the electronic edition of the Astrophysical Journal. A portion is shown here for 
guidance regarding its form and content. 

b The electronic edition of this table contains the continuum fluxes per LH beam area averaged over the wavelength ranges 
10.0 - 10.48 fim and 13.5 - 14.3 fim and the continuum fluxes and errors per LH beam area observed at 15.56, 18.71, 22.93, 
28.22, and 33.48 fim. 
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Fig. I. — Images of the GC region with the observed positions indicated by square boxes. 
Other positions of interest are labeled, (a) Radio continuum (log log scale) imaged at 21 
cm by Yusef-Zadeh & Morris (1987a) with 11" resolution. The Radio Arc and part of the 
Sgr A region are nonthermal; the rest of the radio emission is thermal bremsstrahlung. A 
few positions are labeled with numbers (see Table 1). Also indicated is the size (11" by 22") 
and orientation of the LH slit. Including the small maps, the actual size of the observed 
positions is 11" by 30" for LH and 13" by 14" for SH. (b) MSX Band E (21 fim) image (log 
scale) from Price et al. (2001). 
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Fig. 2. — Three typical spectra. Lines representative of each type of location are marked 
(almost all lines are detected at each position). The broad features at 11.0, 11.2, 12.0, 12.6, 
13.5, 14.2, and 16.3 - 17.5 fim are due to polycyclic aromatic hydrocarbons (PAHs). (a) 
The center of the Bubble (Position 18). (b) GC H II regions (Position 34 in the Arched 
Filaments). The spikes on either side of the [Ne II] 12.8 jum line are ghosts, (c) The diffuse 
ISM at the extreme distance from the Galactic Plane (Position 38). 
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Position 38 36 31 26 21 16 11 6 1 



in i 1 1 1 i i mi i i i 1 1 i i i i in i i 1 1 1 i i ii 




0.3 0.2 0.1 0.0 -0.1 -0.2 -0.3 -0.4 
Galactic Latitude (degrees) 



Fig. 3. — Normalized fluxes for all positions as a function of Galactic latitude, (a) The H2 
S(l) 17.0 /iin line, the H II region and PDR lines of [Si II] and [Fe II] and the high excitation 
[O IV] line are plotted. The normalizing fluxes (SH fluxes scaled to the LH beam = 248 
square arcsec) are H 2 S(l): 1.45 x 1(T 15 W m~ 2 , [Si II] 34.8 /an: 2.33 x 1(T 14 W m" 2 , [Fe II] 
26.0 /mi: 6.93 x 10~ 16 W m" 2 , and [O IV] 25.9 /mi: 1.49 x 10~ 16 W m" 2 . (b) H II region 
lines of [Ne II], [Ne III], and [Fe III] are plotted, along with the 28 /im continuum from the 
PDR-heated molecular clouds. The normalizing fluxes are 28 /an continuum: 97 Jy, [Ne II] 
12.8 /im: 5.09 x 10~ 14 W m" 2 , [Ne III] 15.6 /an: 5.65 x 10~ 15 W m" 2 , and [Fe III] 22.9 
/on: 1.64 x 10~ 15 W m -2 . The [S III] lines are not plotted since they have very similar 
distributions to the [Ne II] line. 
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Position 38 36 31 26 21 16 11 6 1 
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Fig. 4. — Radial velocities Vlsr as a function of Galactic latitude. A few positions are labeled 
with numbers. The velocities of the ionized lines are normalized to —15 km s _1 at Position 
31 in the Arched Filaments, based on radio recombination line measurements there (Lang et 
al. 2001). The dotted line is the average of the [Si II], [Ne II], and [S III] velocities at each 
position (Table 1). The average H 2 S(l) line velocity is normalized to the average of the CO 
J = 4 - 3 velocities measured by Martin et al. (2004) and the solid line is these velocities 
at each position. For most positions, the measured uncertainties for the line centers are < 2 
km s _1 for the brightest lines: [Ne II], [S III], and [Si II], and < 10 km s _1 for [Ne III] and 
H2 S(l). The dispersions in the difference of the measured velocities for pairs of bright lines 
are ~ 10 — 20 km s _1 ; this probably represents the actual measurement uncertainty for the 
brightest lines. 
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Position 38 36 31 26 21 16 11 6 1 
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Fig. 5. — Electron density (N e ) estimated from the [S III] 18.7/33.5 /iin line flux ratio 
plotted as a function of Galactic latitude. The Galactic latitudes of the Arches Cluster and 
the Bubble are marked. Most of the error bars are due to uncertainties in the extinction. 
For 8 of the 10 positions with N e = 10 cm" 3 (Table 1), the [S III] 18.7/33.5 fxm line flux 
ratio computed with the original estimated extinction was lower than the theoretical lower 
limit for T e = 6000 K; for these positions the extinction was increased until a density of 10 
cm -3 could be obtained. 
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Fig. 6. — Ionic abundances ratios plotted as a function of Galactic latitude, (a) The (Ne" 
+ Ne ++ )/H+ ratio, (b) The (S++ + S +3 )/H+ ratio, (c) The Si+/(Ne+ + Ne ++ ) ratio. 
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Fig. 7. — Ratios indicative of gas excitation plotted as a function of Galactic latitude, (a) 
The Ne ++ /Ne + ratio, (b) The S +3 /S ++ ratio. This ratio has large error bars because the 
[S IV] line suffers so much extinction. 
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Fig. 8. — The [Si II] 34.8 /im/[S III] 33.5 fim line flux ratio plotted as a function of Galactic 
latitude. 
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Fig. 9. — Ionization fractions and continuum flux ratios affected by dilution of the radiation 
field, plotted as a function of distance from the Quintuplet Cluster on the plane of the 
sky. The measured ratios are for Positions 4 - 26, excluding Position 13, which may be a 
separate compact H II region. In addition, ratios computed from models calculated using 
CLOUDY (Ferland et al. 1998) are plotted as asterisks connected by solid lines (see text), (a) 
Ne ++ /(Ne + + Ne ++ ). (b) S +3 /(S ++ + S +3 ). (c) Ratio of the extinction-corrected continuum 
fluxes at 18.7 and 33.5 /xm. 
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Fig. 10. — Normalized spectrum segments (continuum subtracted) showing weak but impor- 
tant lines (see Fig. 2 for typical strong lines). The spectra are plotted in order of position 
number (Table 1) from bottom to top and are labeled every 5 positions, (a) The 17.8 - 18.0 
spectral region has the [P III] 17.89 /xm and [Fe II] 17.94 /im H II region lines, (b) The 
25.77 - 26.09 spectral region has the high excitation [O IV] 25.89 /im line and the [Fe II] 
25.99 /im PDR line. 
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Fig. 11. — The +3 /(Ne + + Ne ++ ) ratio plotted as a function of Galactic latitude. 
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Fig. 12. — The +3 /0 ratio plotted versus the Ne ++ /(Ne + + Ne ++ ) ratio from various 
H II region models (see text). An estimate of the observed +3 /0 ratio is also plotted as 
red crosses with error bars, where the estimate is calculated by multiplying the measured 
+3 /(Ne + + Ne ++ ) ratio by the Ne/0= 0.25 ratio found for the Orion Nebula (Simpson et 
al. 2004). The six positions with diamonds overplotted are Positions 14 - 19 in the Bubble. 
The open circles are the CLOUDY models with the Sternberg et al. (2003) atmospheres 
for T c g = 50 kK (no model with a cooler ionizing SED produced enough +3 to appear on 
this plot. The asterisks are CLOUDY models ionized by multiple component black bodies; 
the values of T eS used were 30, 35, and 40 kK for the cool components and 10 5 K (white 
dwarfs) or 10 6 K (X-rays) for the hot components. The triangles are models with ionizing 
SEDs from Sternberg et al. (2003) (T e ff = 36 kK) for photon energies E < 4.013 ryd and 
increased by arbitrary amounts for E > 4.013 ryd (see text). The models represented by 
the four triangles connected by lines all have SEDs with emission in the He ++ continuum 
(E > 4.013 ryd) of a factor of 10 5 larger than the continuum at 4.0 ryd. The four models 
have Ri = 2, 5, 10, and 15 pc (top to bottom). Models with completely unphysical SEDs can 
produce ratios approaching the observed ratios, but these models do not have the decrease 
in excitation with distance from the ionizing stars as is observed. 
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Fig. 13. — The (Fe + + Fe ++ )/(Ne + + Ne ++ ) ratio plotted as a function of Galactic latitude. 
The red crosses with larger error bars use the [Fe II] 17.94 jum line to measure the Fe + 
abundance and the blue open circles with the tiny error bars use the [Fe II] 25.99 /zm line. 
The 25.99 /im line is additionally formed in PDRs at much lower T e than the 6000 K used 
here, whereas the 17.94 /im line is exclusively formed in H II regions like the [Fe III], [Ne II], 
and [Ne III] lines. Thus, the crosses are more likely to represent the abundance ratios of the 
H II region gas. 
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Fig. 14. — The (Fe + + Fe ++ )/Fe ratio plotted versus the Ne ++ /Ne ratio for various H II 
region models (see text). Such models can be used to estimate the correction for higher 
ionization stages of Fe in a determination of the abundance of gas-phase iron. The NEBULA 
models are for H II regions with A^yc = 10 49 s -1 and Ri = 0, whereas the CLOUDY models 
have N Lyc ranging from 10 49 to 10 51 s -1 and Ri ranging from 0.01 to 10 pc. The heavy solid 
line and open circles are the CLOUDY models with Sternberg et al. (2003) supergiant SEDs, 
T eff = 32-37 kK, N hyc = 10 51 s _1 , and Ri = 10 pc (Log U = -2.05). These are the models 
most likely to represent the regions of the Bubble ionized by the Quintuplet Cluster and the 
regions of the Arched Filaments ionized by the Arches Cluster. 



1.0 
0.8 



CO 

S 0.6 



S 0.4 

CM 

0.2 

0.0 
1.5 



CO 

1.0h 

CO 
CM 

ll 0.5 



0.0 







+ 




X 



X 



X 



X 



X 



X 

X' + - 



x 



x 
x 



x 



X 



.* • X 

X 



-x 

X x 



XX 



XX x, 



x^ 



X 



T 



+ + x 



Xh 



X 



X 



4, +X 



X 



X 



2 

X 9.6 



Fig. 15. — Molecular hydrogen ratios plotted versus optical depth at 9.6 /iin, rg.6, where 
rg 6 is determined by fitting the ratio of continuum fluxes at 10 and 14 /iin with a template 
from the Orion Nebula. The crosses with error bars are the flux ratios corrected with the 
extinction-versus-wavelength relation of Chiar & Tielens (2006) for the GC; the asterisks 
are the same data corrected using the extinction law of Draine (2003a, b). The solid and 
dotted lines represent least squares fits through the two sets of extinction corrected data, 
respectively, (a) The H 2 S(0)/S(2) ratio. The solid line has a slope of -0.06 ± 0.02 and the 
dotted line has a slope of -0.12 ± 0.02. (b) The H 2 S(2)/S(l) ratio. The solid line has a 
slope of -0.001 ± 0.003 and the dotted line has a slope of 0.195 ± 0.003. 
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Fig. 16.— The H 2 S(0)/S(1) flux ratio plotted versus the H 2 S(2)/S(l) flux ratio. The crosses 
with error bars are the ratios corrected for extinction using the extinction law of Chiar & 
Tielens (2006) for the GC; the asterisks are the same data corrected using the extinction law 
of Draine (2003a, b). The solid line connects the ratios of the theoretical fluxes of these lines 
in LTE; the diamonds mark the gas temperatures, T, along the line. The dashed line shows 
the locus of ratios of the sum of a cool component with T = 125 K and a hot component with 
T = 1000 K and the dotted line shows the locus of ratios of the sum of a cool component 
with T = 100 K and a hot component with T = 450 K. The squares mark where the ratio 
of the column density of the cool component to the column density of the hot component is 
1000, 100, 10, and 1. 
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Fig. 17.— The H 2 S(0)/S(2) flux ratio plotted versus the H 2 S(2)/S(l) flux ratio. See Fig. 16 
for a description of the lines and labels. 



